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Introduction 


During the summer of 1972 four tests were conducted to provide flight conditions 
for qualifying the Vikiimg Decelerator System in a simulated Mars environment. 
Because of uncertainties in the Martian atmosphere, it was necessary to test 
the decelerator at supersonic, transonic and subsonic conditions. A Balloon 
Launched Decelerator Test (BLDT) vehicle which has an external shape similar 
to the actual Mars Vikimg Lander Capsule was used so that the decelerator 
would be deployed in the wake of a blunt body. The BLDT vehicle and inflated 
decelerator are shown im Figure 1. A detailed account of the flights is 
given in References (1), (2), (3) and (4). 

As part of the post-flight analysis an effort was made to simulate the BLDT 
vehicle flights from dhe time they were dropped from the balloon, through 
decelerator deployment, nxntil stable decelerator conditions were reached. 

The purpose of this paper is to outline the procedure used to simulate these 
flights using the Statistical Trajectory Estimation Program (STEP). Using pri- 
marily ground-based position radar and vehicle onboard rate gyro and accelerometer 
data, the STEP produces a minimum variance solution of the vehicle trajectory 
(position and velocity) and calculates vehicle attitude histories (Euler angles). 
Using film from cameras in the vehicle along with the computer program described 
in Reference (5), attitude histories for portions of the flight before and after 
decelerator deployment were calculated independent of the STEP simulation. With 
the assumption that the vehicle motions derived from camera data are accurate, 
a comparison reveals that STEP was able to simulate vehicle motions for all 
flights both before and after decelerator deployment. 


Tools 


The primary tool for simulation of the BLDT vehicle flights was the Statistical 
Trajectory Estimation Program (STEP); more specifically, STEP 2. A formulation 
report. Reference (6), presents the equations, theory, and numerical techniques 
used in the STEP. A utilization report. Reference (7), presents the computer 
program listing, input and output information, and general user information. 

These reports should be obtained and reviewed before attempting to run the 
STEP. This paper does mot attempt to replace these reports, but is a comple- 
mentary work presenting program utilization for a specific task. 

No basic modifications were made to the STEP for the BLDT simulations. Output 
changes were made so that an output tape could be generated in a format 
compatable with an existing plotting program. A total angle of attack calculation 
was also added during the output phase of the program. 

A complete listing and flow chart of the STEP is given in Reference (7). 

A listing of the control cards necessary to retrieve the STEP, which is on a data 
cell at Langley, and the changes which were made to the program are presented 
in the Appendix. The control cards as listed enable running STEP, program 
TAPES, and a plot program BRUTUS, as one input job. Program TAPES separates 
the STEP output tape into two tapes, one containing the forward integration, 
the other the backward smoothing pass. This program properly formats the data 
for the plotting program. BRUTUS, the plotting program, plots designated 
parameters on the CALCOMP plotter. Listings of these two routines along with 
the necessary information for using BRUTUS are given in the Appendix. The 
listings present the STEP, TAPES, and BRUTUS as they were input for the AV-4 
example simulation. 
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A smoothing program, SNK30TH, was utilized to smooth telemetry data before it was 
used in the STEP. This program uses a least square polynomial curve fit with 
a walking second degree mid-point filter. It incorporates the Langley computer 
facility subroutine LSQP0L which is documented in Reference (8) . 


Implementation 

The transonic and supersonic flights were simulated in two parts. The first 
part is from vehicle drop until the decelerator is deployed (i.e., mortar fire). 
The second starts at mortar fire and simulates the vehicle motions with the 
decelerator deployed. The simulations were divided at mortar fire primarily 
because the rates and accelerations after mortar fire have greater amplitudes 
and higher frequencies than those prior to mortar fire. A secondary reason 
is that inaccuracies in telemetry data or initial conditions are partially 
additive in STEP because of the integration and, therefore, errors may 
increase with time. : 

The basic simulation process is the same for before or after mortar fire. 

Any exceptions are discussed as they arise. The actual simulation process begins 
by accumulating all necessary input data. First, a binary tape containing time 
histories of slant range, azimuth and elevation angle from the primary radar 
along with corresponding standard deviations for these parameters is prepared in 
the manner described in the utilization report. This is called a FIT tape. 

Little difficulty was encountered in preparing the FIT tapes since the radar 
data is not conditioned before being used in STEP. Time histories of slant 
range, azimuth and elevation angle are shown in Figures (2), (3), and (4). 

A second binary tape, the PQR tape, containing the vehicle telemetry rate gyro 
and accelerometer data must also be prepared. This data must be conditioned 
and. a discussion explaining how this can be accomplished is given in the appendix 
of Reference (7). In this instance, the wild points were removed by a difference 
technique similar to that in Reference (7) or, in areas of high non-linearity, 
they were visually recognized and replaced individually. The data were then 
smoothed using the smoothing program already discussed. In some cases, especially 
those simulations from drop up to mortar fire, as many as 29 points of data 
were used in the smoothing filter and as many as three smoothing passes were made. 
Some cases dealing with data after mortar fire were not smoothed because of the 
rapid oscillation of the data. In these instances only the wild points were 
removed. It should be remembered the primary objective of the data conditioning 
is to obtain a time history for each parameter that most accurately reproduces 
the rates and accelerations actually encountered in flight. The order in which 
the PQR data are written on tape is given in Reference (7). Time histories of 
the pitch, yaw and roll rates and X, Y, Z accelerations for the period from drop 
to mortar fire are shown in Figures (5) through (10). Figures (11) through (lb) 
show the rate and acceleration inputs for the period after mortar fire. 

The remaining data are submitted on punched cards. The formats along with an 
explanation of the parameters are given in the utilization report. A sample 
listing of the card input is shown in the Appendix. A brief discussion of each 
input card as jt pertains to this project follows. 

I ! 

The first ; card 1 is a comment card used to label the run. The second card 
contains 10 control parameters. The first parameter, NPC (1) determines the 
mode of operation. The filtering mode was used for this work because it 
was designed for post flight simulation. The second parameter NPC (2) indicates 
whether English or metric units are to be used. NPC (3) determines the form of 
the input and output variables. NPC (3) was set equal to one so that the initial 
attitude as well as the vehicle attitude histories would be in Euler angles. 
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NPC (4) determines whether the reference trajectory as defined by radar will be 
updated or not during the forward integration. Both methods were tried, but 
only the updated mode would give satisfactory results. NPC (5) was set equal 
to zero so the radar data would be printed, NPC (6) is also a print-control 
to indicate whether or not the covariance and correlation matrices are to be 
printed. NPC (7) indicates the number of iterations. As many as three iterations 
were attempted, but the best results were obtained on the first iteration. 

NPC (8) was set equal to 3 to smooth the state variables and covariances and 
calculate residuals and loss function. The residuals were used to determine the 
difference between the original and final STEP radar data. NPC (9) set to zero 
processes radar slant range, azimuth and elevation angle data by the vector 
process which is most efficient. NPC (10) was set equal to zero to prevent a 
STATE tape from being generated. Instead, a tape of output data computable with 
the plotting program was -written in the output portion of the program. Most 
control options were exercised to determine exactly what effect their operation 
had on the STEP. Those given in the sample proved to be the best combination for 
the BLDT simulations. 

The next nine cards sett the initial state parameter values and their standard 
deviations. For the simulations initiated at drop, values for velocity, 
altitude, flight path angle, azimuth, latitude and longitude were obtained from 
radar data. Of the three Euler angles, Psi (vehicle azimuth) was established 
using data from the magnetometer on the load bar. Theta (vehicle pitch attitude) 
and Phi (vehicle roll attitude) were obtained from preflight measurements. 

For those simulations initiated at mortar fire, the initial state variables were 
derived from the final conditions of the best STEP simulations which were 
terminated at mortar fire. The standard deviations of the state variables were 
obtained either from the accuracy of the radar data or, in the case of the Euler 
angles, by engineering judgement. The varying of the initial conditions and 
their standard deviations influenced the trajectory simulation. Those shown in 
the sample are typical of what were necessary to obtain a good simulation. 

The next twelve cards of the sample input are examples of the biases and 
scale factors of the rates and accelerations. Their model parameter number 
and the proper units are shorn in the sample listing as well as in the utiliza- 
tion report. 

The next three cards identify the latitude, longitude, and altitude of the 
radar site. 

The next card specifies initial time, printout frequency, computing interval, 
and final time. The computing interval should be small enough to accurately 
define the rate and acceleration time histories. In the simulations before 
mortar fire a computing interval of .1 second was used. For simulations 
after mortar fire the computing interval was .02 second. 

The next card specifies geophysical and gravitational data. In all the simula- 
tions the polar and equatorial radii were equal to the local earth radius at 
White Sands. This was done to bypass a suspected problem with the oblate earth 
subrout inb . 


5 



The next to last input card identifies the type of tracking data; in this case 
radar. It specifies tie initial and final time to process the radar data, and 
the processing interval!. It also specifies processing of all three radar 
components and an incnanental standard deviation for each component. These 
standard deviations are added to those on the FIT tape to produce a total 
standard deviation. TMs completes all card input data. The last card has only 
the number 20 beginning in Column 1. 

The program should now be ready to run. The primary objective is to simulate 
the BLDT flights by matching vehicle attitude time histories derived from 
on-board camera data vtele still comparing favorably with the radar parameters. 
After the initial STEP run is made, using nominal or default biases and scale 
factors, an attempt is made to match the roll attitude time history obtained 
from the STEP with thatt from camera data. This was accomplished by adjusting 
the scale factor or bias on the roll rate data or, in the simulations after 
mortar fire, by adjusting the initial roll angle. Phi. As in the example of 
the drop to mortar fir® portion, the roll attitude at mortar fire was approximately 
75 degrees larger than it should have been when the default scale factor of 
1.0 was used. By looking at the roll rate telemetry data, the average spin 
rate during the approximately 30 seconds the vehicle is spinning is 150 degrees/ 
second. It can then be calculated that a scale factor of .9835 will decrease 
the total roll by 75 degrees. The program is rerun with the new input and 
compared again. This {process continues until the camera roll attitude data 
compare favorably with the backward smoothing simulation of the STEP. Next, 
the pitch and yaw attitude histories are compared to determine if a bias or 
scale factor adjustment is necessary. As these iterations proceed, make sure 
the roll attitude does not change appreciably due to changes in pitch or 
yaw. Furthermore, as the attitude histories begin to compare favorably, the 
radar parameters derived from the STEP (i.e., altitude, velocity, flight path 
angle, and azimuth) should be checked against the actual radar data. 

Discrepancies in the radar data may also be adjusted with a bias or scale factor 
or by a change in initial conditions. Changes in biases, scale factors, or 
initial conditions should be within the accuracy of the telemetry or radar data. 

Any large changes that seem necessary would indicate a problem on either the 
FIT or PQR tape or possibly an incorrect initial condition. The iteration 
process continues until the STEP . attitude histories compare favorably with those 
derived from camera data and the STEP radar data checks with actual radar data. 

Figures (17) through (19) show the difference between the STEP attitude histories 
and camera data from drop until mortar fire. The comparisons of flight path 
angle and azimuth with radar are also shown in Figures (17) and (18). Comparisons 
of the STEP velocity aasd altitude histories with radar data are shown in Figures 
(20) and (21). Figures (22) through (26) present similar comparisons for the 
period after mortar fixe. 

These comparisons reveal that the STEP was able to simulate the BLDT vehicle 
dynamics and trajectory parameters successfully. It should be remembered that 
an iterative process was used to obtain these results and the steps described 
above are for a specific simulation. For simulation of other vehicles and 
trajectories the order or possibly the necessity of some steps will probably 
be different. This example does, however, provide insight on how the iterative 
process was able to produce good results for this particular problem. 
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Figure 1 BLDT Vehicle and Decelerates 
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Figure 6. Telemetered Yaw Rate from Drop to Mortar Fire 
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Figure 7. Telemetered Roll Rate from Drop to Mortar Fire 
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Figure 8. Telemetered X-Axis Acceleration from Drop to Mortar Fire 
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Figure 9. Telemetered Y-Axis Acceleration from Drop to Mortar Fire 
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Figure 10. Telemetered Z-Axis Acceleration from Droo to Mortar Fire 
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Figure 14. Telemetered X-Axis Acceleration after Mortar Fire 
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Figure 17. Vehicle Attitude and Flight Path Angle from Drop to Mortar Fire 
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Figure 18. Vehicle Heading and Azimuth from Drop to Mortar Fire 
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Velocity from Drop to Mortar Fire 
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Fisure 24. Vehicle Roll Attitude after Mortar Fire 
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Figure 26. Altitude after Mortar Fire 
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SAMPLE COMPUTER INPUT 


i 


JOB , 1 4 0600 , 1 ooono , 08000 « E 1 703 RYJ107 100349 3IN NO* 67 

USER. DOE, J 971561868 67200 NAS1-10900 LTV 

L 1 NECNT < 10000 ) 

REWIND (SCEILF ) 

FETCH ( E 1 71 9 * SPRA06 » SOURCE ) 

NORFL* 

Run ( s * , »sceile, trash > 

REQUEST , TAPE 1 , HY • 

REQUEST , TAPE2 , HY. 

REQUEST , T APE9 * HY • 

REWIND (TAPEl > 

REV! I NO (TAPE2 ) 

RFV'IND(TAPF9) 

sftcore • 

LOO. 

DROPFIl (TAPE l ) 

0R0PF1L (TAPE2 ) 

Rijn (S , , , , .billy, , , > 

REQUEST, TAPEl 0»HY. 

REQUEST. TAPEl 1 , HY * 

R^W I NO ( T A PE 9 ) 

REWIND (TAPE 10 ) 

REWIND (TAPF1 1 > 

setinof. 
billy, 

DROPF IL (TAPE9 ) 

DROPEIL (TAPEl O > 

REWIND (TAPEl 1 ) 

RUNCS. , , , ,M IKE, , , > 

SETINDE* 

MIKE. 

DROPF I L ( TAPE 1 1 ) 

RE W I ND ( CALTPE ) 

REQUEST, TAPFa,Hl ,X. CALTP , R I M . EWS, STEP SIM 

REWIND(TAPF8 ) 

COPYRF (CALTPE. TAPE8 ) 

UNLOAD (TAPES) 

EXIT. 

REWIND (CALTPE ) 

REQUEST, TAPES, HI ,X. CALTP , R I M , EWS , STEP SIM 

REWIND (TAPES ) 

COPYBF (CALTPE, TAPE8) 

UNLOAD (TAPES > 

DROPF IL (TAPEl 1 > 


71 4007, ROL , TEC, RADAR-STEP FOR AV-4 
4 I 305A ,ROL ,EWS, AV-4 STEP? INPUT PQR 
4060 17,RIL,EWS, STEP PLOT 


404033 »RIL »EWS , STEPPLOT FORWD AV-4 
404039, RIL, EWS, STEPPLOT BCKWO AV-4 



DPOPF 1 L ( TAPE 2 5 
DROPF !L ( TAPF9 > 

DROPF IL (TAPF1 0 ) 

DROPF !L (TAPF1.1 ) 

< 789 CARO) 

CUTOUT 2 

2 TAPE1 ,TAPE2»TAPE3,TAPE4*TAPE7,TAPE9> 

CUTOUT 50 51 

CUTOUT 61 62 

CUTOUT 78 

44 CALL OUTPT 
CUTOUT 701 

IF (EOF, IN ) 2000 ,2001 

2001 WRITE (OUT, 1001) KDUM , < B ( I ) , I = 1 » 8 > 

PUTIN 864 

2000 STOP 
CUTOUT 1005 

SUBROUTINE OUTPT 


PUTIN 

1 843 





I F < ST ( 1 ) + CT < 1 ) 

• EO. 

0* 

) CT ( 1 ) = ,0000 

PUTIN 

1847 





I F ( ST ( 2 ) +CT ( 2 > 

• 

c 

UJ 

« 

0. ) 

CT (2 )®. 00001 

PUTIN 

1869 





IF <ST(4 )+CT(4 ) 

• EQ* 

0. ) 

CT (4 )=. 00001 

PUTIN 

1 873 





I F ( ST ( 5 ) +CT ( 5 ) 

• EQ* 

0. ) 

CT(5)=*00001 

PUTIN 

1 880 





I F (ST ( 3 ) +CT ( 3 > 

• EQ* 

0* ) 

CTO )=. 00001 

PUTIN 

1 882 





IF(SP < 1 )+VO( 1 ) 

• EG. 

0* ) 

VO( 1 ) = *00001 

CUTOUT 

1896 





V0(t7)=V0(17)/57.3 
V0(18)=V0< 181/57.3 

BILLE=ATAN(SQRT(TAN(VO( 17) )#*2+TAN< VO ( 18) )**2) ) 

BILLE=BILLE*57«3 
V0< 1 7)=V0( 17)*57.3 
VO ( 18 ) = VO ( 18 )*57*3 

• WRTTF(OUT, 1008) ( XP ( T ) « I = 1 . 3 ) , ( SP < I ) * I = 1 , 2 ) ,B I LLE 

PUTIN 1899 

ppppp=PA (1 )*57*296 
OOOQO=PA (2 )*57.296 
RRRRR=PA ( 3 ) *57 • 296 

WR ITE(9 )TO, VO < 1 ) , VO (2 ) ,VO ( 3 ) ,VO(4 ) , VO ( 17) , VO( 18 ) , PPPPP, QQQOO ,RRRRR 
1 , AB < 1 ) , A8 ( 2 ) * AB ( 3 )', VO (10), VO (11), VO ( 12 ) , B I LLE 



CUTOUT 2839 

24 CALL OUTPT 
(789 CARD) 

1 Ay- 4 f.TEP SIMULATION BEFORE MORTAR FIRE 


2 


0 

1 1 o 

c 

o 

o 

3 0 

0 

3 

l 

1 

100*1 

2. 

- 

INITIAL VELOCITY 

3 

1 

2 

! .4 

.020 


INITIAL FLIGHT PATH ANGLE 

3 

1 

3 

-126. 

.02 0 


INITIAL AZIMUTH ANGLE 

3 

-1 

4 

I 20762. 

50. 


INITIAL ALTITUDE 

3 

1 

5 

3 3.28710 

.001 


INITIAL LATITUDE 

3 

1 

6 

-106*23220 

.001 


initial LONGITUDE 

~3 

~1 

7 

-137. 

.0001 


VEHICLE ROLL ATTITUDE* PS I 

3 

1 

8 

55. 

.0001 


VEHICLE PITCH ATTITUDE* THETA 

3 

1 

9 

0. 

.01 


VEHICLE YAW ATTITUDE* PHI 

3 

5 

36 

.9835 



ROLL SCALE FACTOR 

3 

5 

4° 

1 • 



PITCH SCALE FACTOR 

3 

5 

44 

1 • 



YAW SCALE FACTOR 

3 

5 

45 

-.01745 



ROLL BIAS 

3 

5 

46 

-.010646 



PITCH BIAS 

3 

5 

47 

.. -.06108 



YAW BIAS 

3 

5 

61 

1 • 

• ■■ ■ 


X-AXIS ACCELERATION SCALE FACTOR 

3 

5 

65 

1 . 



Y-AXIS ACCELERATION SCALE FACTOR 

3 

5 

69 

1 • 



Z-AXIS ACCELERATION SCALE FACTOR 

3 

5 

70 

0. 



X-AXIS ACCELERATION BIAS 

3 

c 

71 

0. 



Y-AXIS ACCELERATION BIAS 

3 

5 

72 

0 . 



Z-AXIS ACCELERATION BIAS 

3 

5 

88 

32.90139333 


- V V ■-* .< <■. w 

• P-123- RADAR" SITE- LATITUDE 

3 

5 

B9 

-106.098614 



R— 123 RADAR SITE LONGITUDE 

3 

5 

90 

4149.53 



R-123 RADAR SITE ALTITUDE 

4 



.2 

• 2 

.1 

o 

• 

c 

<t 

5 

2 

20909 I 58*8 

20909158.8 

.00007292116 

7 

1 


0.20 40 

.4 0 

10 1 1 

1 6.666 .00001333 .00001333 

20 







(789 

CARD) 






PROGRAM TAPES < I NPUT « OUTPUT « T APE5 = I NPUT , T APE 6 = 0 UTPUT .TAPE9 , TAPE 1 0 ♦ T 
1 APF1 ) 5 

- DIMENSION A(25)«AA(25) 

199 FORMAT ( 1 28H1 TIME VELOCITY FPA AZIMUTH ALTITUDE BETA ALP 

1 HA ROLL PITCH YAW AX AY AZ PS I THETA PHI 

2 ETA ) 

299 FORMAT ( 1 28H SEC FT/SEC DEG DEG FT DEG DE 

IG DEG/SECDEG/SECDEG/SEC G.S G.S G.S DEG DEG DEG 

- - 2 DEG / ) 

DO 10 1=1,136 
READ<9 ) < A (L ) «L=1 , 1 7) 

A ( 1 1 ) = A ( 1 1 >/32. 174 
A< 1 2 )=A ( 12 )/32. 174 
A ( 1 3 >=A ( 1 3 J/32. 174 
WRITE U 0 > (A <L ) ,L=1 ♦ 1 7) 

10 CONTINUE 

END FILE 10 

DO 20 1=1 , 135 

READ ( 9 M A A <L) *L = 1 *17) 

20 CONTINUE 

DO 30 1=1,136 

READ (9 ) (AA (L > .L=I ♦ I 7 > 

AA U I >=AA < 1 1 )/32, 174 
-e, AA ( 12 )=AA ( 12 )/32, 174 

AA ( 13 >=AA ( 13 1/32.174 
WRITE < 1 I ) (AA <L ) *L=1 * 17 ) 

BACKSPACE 9 
BACKSPACE 9 
30 CONTINUE 

END FILE 1 1 

100 FORMAT ( 1F7.2»1F10.2»2F7.2,1F12»2*12F7.2> 

STOP 

END 

<789 CARD) 



1 

2 — 

24 

25 


28 

27 

29 

30 

4 

45 

5 


6 

61 
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PROGRAM BRUTUS ( I NPUT * OUTPUT * TAPE 5 = I NPUT « T APE 6 = 0 UTPUT . TAPE 1 1 . TAPES) 

DIMENSION 1 J ( 17 ) . A ( 1 7 « 300 ) »X (300 > . Y (300 5 

0 I MENS ION I PH { 1 00 > « BCOX ( 8 ) . BCDY < 8 ) 

CALL CALCOMP 

CALL LEROY 
IT =0 

CALL CALPLT ( 0 • » 0 . *— 3 ) 

FORMAT ( 215 ) 

READ ( 5 ♦ 2 )NI . ( I J (N ) .N = 1 .N I) 

— FORMAT (1615) 

READ ( 5. 2 INCASES* < IPH(M).M= 1 .NCASES) 

1 M= 1 

I PHASE= I PH ( IM) 

J I =0 

DO 27 1=1.1 0000 

READ ( 1 1 ) < A (N« I > .N=l «NI > 

IF(EOF, 1 1 ) 28.27 
JI=JI +1 
1=1-! 

IF ( JI .GE. IPHASE ) GO TO 29 
CONTINUE 
I T = I 
IM= 1 M +1 

READ(5,2)NOW.NGO 
I F ( NOW . £Q . 3 ) GO TO 25 

- <• r ■ a,» **■' t ‘ V, 4 ;• a, *£**:*:_.* . ;S:ii ■ v- ‘ ■*. a v .. y.t ■, 

IF (N0W.EQ.2 )G0 TO 65 

READ ( 5 . 45 ) IL IM.BOTLIM.UPLIM.XNEW 

FORMAT ( 1 I 10.3F1 0.5 ) 

READ ( 5 . 5 ) I X.X IN.XM IN.XMAX. XMLTPY.TMA JX 
READ < 5 . 5 ) IY.YIN.YMIN « YMAX. YMLTPY . TMA JY 
FORMAT ( 1 I 5 . I F 5 . 2 . 4 F 10 . 5 ) 

I F ( NOW . EQ .4 ) GO TO 61 
READ (5. 1 ) NX ♦ NY 
MX=NX/1 0.+.9 
MY =NY/ 1 0. + .9 

READ ( 5 . 6 ) (BCDX (N) .N=l .MX) 

READ( 5 . 6 ) (BCDY<N) .N=l .MY) 

FORMAT ( BA 10 ) 

DO 7 N= 1 iNI 

IF ( ILIM.EQ. IJ(N> ) I T I ME=N 
IF( IX.EO.IU(N) )IXN=N 
IF ( IY.EQ. 1 U(N) > IYN=N 
CONTINUE 

ymin=ymin*ymltpy 



YMAX = YMAX*Yf.iLTPY 
XM I N = .;Ml N#XMLTPY 
XMAX=XMAX*XMLTPY 
DX = (XMAX-XMIN)/XIN 
DY= (YMAX-YMIN )/Yl N 
65 K = 0 

DO 8 1=1 » IT 

IFCAC ITIME. 1 I.LT.BOTlIM) GO TO*' 8 
IF(A< I TIME* I )*GT.UPLIM) GO TO 8 
K = K+1 

X (K ) =A ( 1XN, 1 >*XMLTPY 
Y(K)=A(IYNiI)*YMLTPY 
IF (X(K ) .GT.XMAX ) X(K)=XMAX 
IF<Y(K ).GT*YMAX) Y(K)SYMAX 
IF(X(K)*LT.XMIN) X(K)=XMIN 
IF (Y IK ) .LT.YMJN ) Y ( K > = YM IN 
X<K)=(X(K)-XMIN)/DX 
Y(K > = (YCK >-YMIN)/DY 

8 CONTINUE 

I K =K 

IF <N0W.EQ.2.0R.N0W.EQ.4> GO TO 85 
CALL CALPLT (XNEW,0. ,-3) 

TM= I ./TMA JX 

call axes ( 0. *0. .0. .XlN.XMIN* DX , TMA JX « TM . BCDX » • 1 5»-NX) 

TM=1./TMAUY . . v ..,.v 

CALL AXES (0. ,0. ,90. * YIN, YMIN.DY, TMAUY,TM»BCDY* *15»NY) 

CALL CALPLTtXIN.O. ,3 ) 

CALL CALPLT(XIN,YIN,2) 

CALL CALPLT (0. , YIN, 2 ) 

85 CALL CALPLTCXCl >,Y< 1 ),3) 

DO 9 I =2 * IK 

9 CALL CALPLT (X(I)«Y(I)»2) 

92 GO TO (25,97* I 0.30) NGO 

97 REWIND 11 

IT=0 

GO TO 24 

10 CALL CALPLT(0.,0. ,3) 

CALL CALPLT (0..0. ,999) 

STOP 

END 
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INPUT INSTRUCTIONS FOR PLOTTING PROGRAM 


COLUMNS 

PARAMETER 

DEFINITION 

UNITS 

FORMAT 

1-5 

NI 

Number of A-arrays on the data tape (e.g.: 3) 


1615 

-6=10 

IJ (1) 

A-array indicies ordered as on the data tape 



11-15 

10 (2) 

(e.g.: 100 


l 

.etc. 

etc. 

108 





120 



51-55 

IJ (10) 




1-5 

NCASES 

Number of trajectory cases on the data tape 


115 


IPH (1) 

(e.g.: 25) (100) 



6-10 

Number of phases of each of these cases 



11-15 

I PH (2) 

Ordered as on the data tape (e.g.: 4) 



etc. 

etc. 




IPH (15) 

Use as many cards as necessary filling 



76-80 

from the left most columns (1-5) 



1-5 

N0W 

Operation Instruction for Plotting This 


215 



Data Group: 





= 1, Plot this data group on a new grid 
= 2, Plot this data group on the previous grid 
= 3, Skip this data group 
= 4, Special patch plot option 



6-10 

NG0 

Operation Instruction for After This Plot: 




= 1, Continue to next data group 
= 2, Return to first data group 
= 3, Stop 

= 4, Continue with this same data group 


If Now = 
If Now = 
If Now = 

1-10 

1 continue with Cards 4, 5, 6, 7 » and 8 

2 or 3 repeat Card 3 for next data group 
4 read only cards 4, 5 and 6 

• ILIM A-array index (on data tape) which is used 

1110 

11-20 

BOTLIM 

to limit the plot range (e.g.: = 100 is time 
1 imited) 

Lower limit value of ILIMIT (e.g,: 0. seconds) 

1 FI 0.5 

21-30 

UPLIM 

Upper limit value of ILIMIT (e.g.: 40.3 seconds) 

1 FI 0.5 

31-40 

XNEW 

Absissa value of origin for this grid. 

1F10.5 



Ol-V 


CARD 

5 

6 

7 

8 
9 

NOTE: 


COLUMNS 

I- 5 
6-10 

II - 20 
21-30 
31-40 

41-50 

I- 5 
6-10 

II - 20 
21-30 
31-40 

41-50 

1-5 

6-10 


1-80 


1-80 


If NGO = 1 
If NGO = 2 
If NGO = 3 
If NGO = 4 


PARAMETER 

IX 

XIN 

XMIN 

XMAX 

XMLTPY 

TMAJX 

IY 

YIN 

YMIN 

YMAX 

YMLTPY 

TMAJY 

NX 

NY 


BCDX 


BODY 


DEFINITION 

A- array Index for absissa 

Length of absissa axis (e.g.: 9 inches) 

Origin value of absissa (e.g.: o seconds) 

Maximum value of absissa (e.g.: 45 fecondl) 

Absissa multiplying facto? to scale axis > 

_ label between .001 and 999. ( e q • 7 ) 
Absissa numbered label spacing (e.g!?’l. inch) 

A-array index for ordinate 

Length of ordinate axis (e.g.: 8 inches) 

Origin value of ordinate (e.g.: 0 feet) 

Maximum value of ordinate (eq * annnnn * 
Ordinate multiplying falter « scaS fSSt) 

Ordin te* iSSSVH td.g.: .001) 


Number^of digits in the absissa axis label 
Number^of digits in the ordinate axis label 

Absissa axis label identification (include 
SECONDS) f6CtS) (e ' 9,: TIME ™ LAUNCH 


UNITS 

FORMAT 

inches 

115 
1F5.2 
1 FI 0.5 
1F10.5 
1F10.5 

inches 

1 FI 0.5 

inches 

115 
1F5.2 
1F10.5 
1F10.5 
1 FI 0.5 

inches 

1F10.5 


115 


115 


Ordinate axis label identification (include 

scaled effects) (e.g.: ALTITUDE (1000 FEET) ) 

Return lo c!r d 3 for the fim^aL 5 ™^ 00 th * tape 
no more input cards are required 9r ° UP ° n the tape 
return to Card 3 for the same data group 


8A10 


8A10 



